Analytical analysis of ground states on O-tt long Josephson junctions. 
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We investigate analytically a long Josephson O-yr-junction with several and tt facets which are 
comparable to the Josephson penetration length Xj. Such junctions can be fabricated exploiting (a) 
the d-wave order parameter symmetry of cuprate superconductors; (b) the spacial oscillations of the 
order parameter in superconductor-insulator-ferromagnet-superconductor structures with different 
thicknesses of ferromagnetic layer to produce or tt coupling or (c) the structure of the corresponding 
sine-Gordon equations and substituting the phase 7r-discontinuities by the artificial current injectors. 
We investigate analytically the possible ground states in such a system and show that there is a 
critical facet length Uc, which separates the states with half-integer flux quanta (semifluxons) from 
the trivial "flat phase state" without magnetic flux. We analyze different branches of the bifurcation 
diagram, derive a system of transcendental equations which can be effectively solved to flnd the 
crossover distance Oc (bifurcation point) and present the solutions for different number of facets and 
the edge facets length. We show that the edge facets may drastically affect the state of the system. 
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I. INTRODUCTION 

Due to the recent progress in technology it 
is now possible to fabricate different types of 
TT Josephson junctions (JJs): high- Tc tri-crystal 
grain boundary JJsi, YBa2Cu307-Nb zigzag ramp 
JJ^, Superconductor-Ferromagnet-Superconductor 
(SFS)'^'^ or Superconductor-Insulator-Ferromagnet- 
Superconductor (SIFS)^. 7r-junctions are very promising 
elements for Josephson electronics. It was already 
suggested that they can be used in analog^ and digital^'® 
circuits in classical regime and for implementation of 
qubitsS. 

In this paper, we focus on long Josephson junctions 
(LJJ) consisting of several and tt parts (facets). We 
will call such junctions O-tt-LJJs. LJJs consisting of very 
short (and random) and tt facets, which are naturally 
formed in 45° high- Tc grain boundaries, were studied by 
R. Mints and coauthors in a series of works (see Ref.hd 
and refs. therein). We are more interested in facets with 
the length a comparable to the Josephson penetration 
depth Xj like in artificially prepared structures. These 
are the sizes which will be used in potential devices based 
on fractional vortex dynamics, both in classical and in 
quantum oneaiii^. 

It was foundi^ii^ that at the point where and tt 
facets join, a new type of non-linear excitation may ap- 
pear. This new non- linear solution of the properly mod- 
ified sine-Gordon equation looks like a vortex and con- 
tains one half of the flux quantum and therefore is called 
"semifluxon" . The semifluxon (SF) is always pinned at 



the joining point between and tt facets. The presence 
of SFs was demonstrated experimentally^'^ by scanning 
SQUID microscopy on YBa2Cu307-Nb ramp zigzag LJJs 
in the long facet limit, i.e., when the length of the facets 
a ^ Xj. SFs were also experimentally observed in the 
tri-crystal grain boundary LJJsiSiiLiS 

In this work we study analytically the ground states 
of O-TT-LJJ with arbitrary number of alternating and 
TT facets and the effect of edge facets on the type of the 
ground state. As it was shown earlier for some particular 
cases, there is a critical facet length Oc which separates 
the domains with the two most natural lowest energy 
configurations: the flat phase state and antiferromagnet- 
ically (AFM) ordered array of semifluxons. 

The joining points between and tt facets we will call a 
phase discontinuity points since, e.g., in YBa2Cu307-Nb 
zigzag LJJs, the Josephson phase (/){x) is 7r-discontinuous 
at these points. In the other types of junctions, e.g., SFS 
and SIFS, the Josephson phase is continuous, but one can 
finally arrive to the same equations making a proper sub- 



stitution of variables: (j){x,t) = fj,{x,t) 
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Follow- 



ing Ref. '14' and regardless of the LJJ type we will denote 
discontinuous phase as (f>{x), while continuous (magnetic) 
component of the phase as n{x). 

In section ^] we introduce the model and represent 
general solution ii{x) for arbitrary distribution of tt- 
discontinuity points and present several examples. In 
section Unl we calculate the crossover distances Oc (cor- 
responding to "AFM ordered SF chain" -"flat phase 
state" transition) for N equidistantly distributed tt- 
discontinuity points with the arbitrary length b of edge 



Typeset by REVT^ 



2 



facets. Finally section ITVl summarizes our results. 



A. General stationary solution for finite length JJ 



II. MODEL AND GENERAL STATIONARY 
SOLUTION 



We consider finite length Josephson junction with TV 
"O-tt" conjunction points (7r-discontinuity points). Let 
the coordinates of these points hexj,j — l,...,N, and 
the coordinates of the two ends of LJJ be xq and xn+i- 
We will write and solve all equations in terms of the mag- 
netic component of the phase /i(a;) which is a continuous 
function^'*. Let ^■' (a;), j = 0, . . . , A'', be the piece of nix), 
inside of j-th interval Xj < x < Xj^^i. Thus in total 
we have + 1 intervals enumerated by j = 0, 1, . . . , A^. 
The /J,^", n = 0, 1, . . . correspond to "0" intervals, and 
^2n+i^ n = 0, 1, . . . correspond to "tt" intervals. While 
we use superscript j to denote a piece of the function 
at j-th interval, we use the subscript j to denote the 
value of the function fi{x) at x = Xj, i.e., fij = f-i{xj), 
j = 0, . . . , N + 1. Since we will look for a continuous 
solution /ij are uniquely defined. 

In these notations the sin-Gordon equation reads^^: 



lit = (-IF sm/i^ 



(1) 



Later on, we will need the above time dependent equation 
for analysis of stability of some solutions. At this moment 
we write only the stationary version of Eq. which has 
the form 



f^xx 



{—ly sin fi^ . 



(2) 



It may be integrated once: 



ifii) =C,~2{-^ycos^^^, j = 0,l,...,iV, (3) 



where Cj are integration constants for j-th interval. We 
look for solutions with arbitrary, but equal boundary 
conditions ^x{xo) — fJ-x{xN+i)- Since is the mag- 
netic fielcU^, we may consider either uniform field or non- 
uniform one. In the first case the boundary conditions 
at both ends should be equal as written above. For non- 
uniform field in addition to having iixixa) 7^ fJ-xixN+i), 
we have to include the hx{x) term^'* in the Eqs. and 
For the sake of simplicity in this paper we will con- 
sider only uniform fields. 

End points xg and xjv+i may be either finite or infinite. 
Infinitely long JJ is considered as limit of finite JJ. First 
we construct solutions for finite length LJJ. 



We require that the derivative is continuous at a; = 
3+'^{xj). j = 1,...,N. Then 



Xj, i.e. fiiixj) 



Co = M°(2;o)^ + 2cos/io, 

J 

= /i°(a;o)^ + 2cos/io + 4^(-l)^cos/i,; 
j = l,...N, 



(4) 



CjV = flxi^N+l) +2(-l) COS^AT+i 

The above system involves two different expressions for 
Cn, which produce the first relation among 11 j: 

N 

2cos/io -t- 4^(-l)' cos/Xi = 2(-l)^cos^Ar+i. (5) 
i=i 

After the integration of © one gets inside of the j-th 
interval: 



Xj =± 



dv 



27rnj +( 

{fj, — ajTT — 27rnj )/2 



y^Cj +2 cos(i^ + (Tjir) 
dv 
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1 



2 • 2 
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a,- 



2T:nj)/2,af\ 
2 

± 



(6) 
(7) 



where Uj are some integers which may be different for 
each particular interval; F(x,to) is the Elliptic Integral 
of the First Kind; the lower limit of integration in the 
first integral is convenient for representation of the final 
result in terms of elliptic functions. To write such limits 
of integration we use a;* ^ x„ in the above equations. 
The integration constants x* can be expressed in terms 

of/i, {3 = {),..., Ny 



X* = Xj — aj F[{l^j — (JjTT — 2nnj)/2, a^], 



(8) 



where the sign in aj is taken from the condition that x 
increases when one goes along the junction. It is altering 
in each extremum of the function /i(a;); the variable a is 
such that (T2n = 1, o'2n+i = 0, 71 = 0, 1 , . . . , iV. Functions 
fi^ (x) can be expressed in terms of elliptic functions from 
the Eq. ®: 



2 am 



(9) 



where am(a:,m) is the Jacobi amplitude. The values fj,j, 
j = 0, . . . ,N+1, are the solutions of the following system 
of TV + 1 equations (j = 0, . . . , iV) 



/ij-l-i — 2Tmj + (TjTT + 2 am 



a. 



(10) 
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and Eq. 

Below we will need extremum values of /i^ (if any) : 



fil-^ — ± arccos 



(11) 



Remember, that the fJ--'{x) can be either the piece of 
monotonic or periodic function depending on Cj, namely: 

\Cj \ < 2 : fi-' (x) piece of periodic function (12) 
|Cj| > 2 : II-' (x) piece of monotonic function (13) 

From the above one can get the following restriction on 
the possible values of the parameters rij-. if ^^(xj) > 0, 
then > tij] if fi-!^{xj) < 0, then rij+i < rij. 

Thus Eq. © together with Eqs. ® and JTHIl repre- 
sents general formulae for all possible ground states in 
O-TT-LJ J of finite length with arbitrary number of discon- 
tinuity points and uniform magnetic field. Each partic- 
ular ground states is characterized by the parameters rij 
and Hj. The later are solutions of the system (fllUl . We 
emphasize, that this system may be either solvable or 
not, which depends on values N and positions of discon- 
tinuity points Xj. 



B. General stationary solution for infinitely long J J 

In this section we give some remarks regarding station- 
ary solutions for infinitely long J J. The main difference 
between finite and infinite length is related to the edge 
facets. Let us consider the limit a;o — *■ — c», xn+i — > oo. 
We admit the following zero field boundary conditions: 



lim — lim /i = 0, 

X — ^±oo X — > — oo 

lim fi = 27771^ + (1 — crjv)7r 

X — >oo 

Equations l@J-l(Sl now have the form: 



Co — Cat — 2, 



(14) 



(15) 



Q = Q_i+4(-l) 



2 + 4^(-l)*cosAt,; j = l,...N. 



i=l 



Eq. (O is reduced to the following one: 



N 



^(-l)VosAi» = 



(16) 



Expression ((TT)l for as well as Eqs. (15)1- 1(117)1 stay the 
same for the inner intervals. But for the edge facets (j = 
0, TV) Eqs. © should be replaced with: 

fi^{x) — 4arctan [Go exp (a; si)] ; (l^a) 
^l^{x) = (l-a-jv)7r 

+ 4 arctan [Gjv exp (a; sjv)] + 27rnjv, (17b) 



where parameters si and sjv may be either 1 or (—1) de- 
pending on boundary conditions at infinities l|14f) . They 
define whether the function is increasing or decreasing. 
Integration constants Go and Gat are defined by the equa- 
tions 



/(a;i) 
^^(a;Ar) 



4 arctan [Go exp {xi si)] ; (18a) 

(1 - (TAr)7r 

4 arctan [Gjv exp (xjv sn)] + 27rnjv(18b) 



C. Examples of particular solutions 

The problem to construct solutions is reduced to solv- 
ing the system (O and lfTU)l . We consider examples of 
two types of solutions: the so-called AFM state "f It |and 
the state ttii for LJJ with equidistant distribution of 
discontinuity points x-m^, zero field boundary conditions 



fJ-xixo) = f^xixN+l) = 

and Uj = for all j. Thus < < 2tt, 



arccos 



(-IPQ 



(19) 



(20) 



In this case expression for aj in the Eqs. (|H))- (|13|I can be 
given in the following form: 



2sgn ^xjxj) 
yG~+2 



j = 1, ... ,7V, ao = - 



= (21) 



Following Ref . 0, we consider equidistant distribution of 
discontinuity points Xj = a{j — = 1, . . . , N, but with 
end points xq = —b and xn+i = a{N — 1) + 6, where b is 
the length of the edge facets which may be different from 
a. Due to the symmetry one has 

For even N, j = 0, . . . , iV/2 : 

= fiN+i-j, (22a) 
For odd iV, j = 0,...,(iV-l)/2: 

= TT - fiN+l-j, M(Ar+l)/2 = 7i'/2. (22b) 



a. AFM state in finite length LJJ. For AFM state, 
fi^ should have extremum value inside its interval, which 
means that /i-' is a piece of periodic function l|12|) , so that 
one needs to provide Cj < 2 for all j. The plots of (f){x), 
fi{x), fixix) for a = 2, 6 = 5 and two different values of 
N {N ^3 and iV = 4) arc shown in Fig. [T] 

b- Ttii state in infinite LJJ . Let us consider in- 
finitely long JJ with iV = 4. For this state should have 
extremum value inside of the middle interval: Gi, G3 > 2; 
Go = G4 = 2, and < G2 < 2. Parameters in Eqs. ifTT)) 
and (|18|) are defined as follows: (Tn = 1, si = 1, S2 = — 1. 
The plots of 4'{x), ^i{x), ^x{x) are shown in Fig. |21 
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FIG. 1: The state with AFM ordered SFs. Graphs of (^(2:), fi{x) and fi^{x) for a = 2, & = 5; (a) iV = 3, /lo = 0.0172, 
fii = 1.2371 and (b) N = 4, fio = 0.016, m = 1.1569, = 1.3775. 
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FIG. 2: The state TTii- Graphs of <t>{x), fi{x) and 

for infinitely long JJ with a = 4 and N = 4; fii — 1.6165, 

fj.2 = 4.6086. 



III. CROSSOVER DISTANCE 



In this section we study the transition between flat 
phase state and AFM ordered SF state in O-tt-LJJ with 
equidistant distribution of 7r-discontinuity points \xj+i — 
Xj\=a,j = l,..., N—1 expressing the length of the edge 
facets in terms of a: b = (3a. For this case < fj, < n. 

In the next subsection IIII Al we study the stabiUty of 
the flat phase solution and derive the system of algebraic 
equations for calculation of Qc- In the subsection llll BI we 
study the existence of AFM solution and derive another 
system of equations defining Oc. 



A. Stability of flat phase state 

We study solutions of the time dependent Eq. which 
have small amplitude oscillations around the flat state ^c- 
As it follows from Eq. jSJ, the value of phase fj,{x) = fic 
in the flat phase state can be either or tt. Note that, for 
odd N the symmetry conditions Ij22b|) result in ^ = 7r/2 
at the middle point x — x-^qji. Since [i(x) = const in the 
flat phase state, fi should be equal to 7r/2, but this is not 
a solution of Eq. (^). Thus, we conclude that for odd N 
the flat phase state cannot be realized, so we formally 
take ttc = 0. Below we consider only even N. 

Let us introduce a new function fl — fi — ^ for 

which equation has one of the forms given below 

AL-Mi = (-l^A, Mc = 0, (23a) 
= (-ir+V, /^c = 7r (23b) 

for all j = 0, ...,7V + 1. 

First, we study only the Eq. (|23a|l corresponding to 
He — 0- To study the stability we look for the solution in 
the following form: 

/IJ^e^V^ (24) 

where E is considered around E = 0, because this is the 
point where the stability of solution changes. Hereafter 
we use — 1 < i? < 1. Then Eq. H23a|l gets the form 

i^L^[E+{-iy]i^', J=0,...,N + 1. (25) 

There are two solutions of the correspondent character- 
istic equation: k = ±fci for even intervals and k = ±ik2 
for odd intervals, where ki = y/1 + E and k2 = \J\ — E 
are both real. The solutions of the Eqs. can be rep- 
resented by the following system 

= ^ocosh(fcij;-f /3o) (26a) 
= A, cosh[A:i {x ~ a{j - 1)) + /3,], j = 2,4,. .(^6b) 
i^^ = Aj cos[/c2(a; - a{j - 1)) -f (3,], j = 1, 3, . . . (26c) 
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The zero field boundary condition /j,x{—f3a) = gives the 
formula for f]^: Po = /3a. Due to the symmetry (I22b|l . 
it is enough to consider only half of the whole LJJ. We 
have for the middle interval fJ.N/2 = which gives 

expressions for /3^/2- — — fcia/2 for even N/2, or 

Pn/2 — — fc2a/2 for odd N/2. Continuity of the functions 
/i imposes the following relations among parameters Aj 
and Pj: 



Aq cosh(/3a) 
Aj cos(fc2a + Pj) 
Aj cosh(fcia + Pj) 



Ai cos(/3i 



Aj+icosh(/3j+i), j 
Aj+i cos(/3j+i), j = 



(27a) 
= 1,3, (2.7b) 
2,4, ..(.27c) 



To provide continuity of , one needs to impose addi- 
tional relations among the parameters Pj: 



ki tanh(/3a) 
/c2 tan(fc2a + Pj) 
ki tanh(fcia + Pj) 



-/C2tan(/3i), (28a) 
-fcitanh(/3,+i), j = l,3,.(2^b) 
-fc2tan(/3,+i), j = 2,4,..(?8c) 



where j < N/2. Equations (|27|l define the amplitudes 
A-j, j > 0, in terms of Aq and /3„, n = 0, . . . , N/2. The 
Eq. 128a(l establishes relation between a and E, while 
Eqs. (|28b|l and l|28c|l define parameters Pj. 

The solution ^ is stable if < 0. We define 
crossover distance Uc the distance for which E ~ and, 
consequently, ki — k2 — 1, i.e., the system of Eqs. (I28a|l - 
^cjl defining Oc is 



tanh(/3ac) 
tan(ac + Pj) 
tanh(ac + Pj) 

Pn/2 



■tan(/3i), 

■tanh(/3j+i), j 1,3, . 
-tan(/3j+i), j-2,4,.. 



(29) 



Note that this system of equations is rather easy to solve 
consequently excluding Pj. At the end one gets a tran- 
scendental equation for Uc [the first equation of if^ ]. 
with given parameter P and function Pi(ac). It is clear, 
that flc = and all Pj = is solution of the system 
(|29|l . Instead, we are interested to find the first non-zero 
solution flc of (|29|l . 

Non-zero solution to this transcendental equation does 
not exists for any p. To derive the existence conditions we 
analyze the behavior of the function tan[— /3i(a)]. From 
Eqs. H29|l it is evident that tan[— /3i(a)] is monotonic func- 
tion of a with uniform concavity. This conclusion is direct 
consequence of the fact that both tan(x) and tanh(a;) in 
the system (|29ll are monotonic with uniform concavity 
inside of the continuity interval < a; < 7r/2. Thus, if 
tan[— /3i(a)] increases in some particular point, then it 
will be increasing inside the whole interval. Let us find 
the behavior of tan(— /3i(a)) for a ^ 0. All but first equa- 
tions of the system (|29f) are easily resolvable for Pj , j > 1 
giving Pj « —a/2. Thus tan[— /3i(a)] is an increasing 
function of a. Since tan(a;) — > oo for x —^ 7r/2, this func- 
tion is concaved up. Since tanh(/3a) is concaved down, 
the two functions may intersect [the first of Eqs. H29|) 
has nonzero solution], only if tanh'(O) > tan'(/3i(0)) (the 







N 


(3 = 1/4 


P = l 


/3 = 2 


P = oo 


2 


2.92771 


1.4639 


1.5670 


tt/2 


4 


1.25461 


1.1772 


1.2989 


1.3063 


6 


0.98327 


1.0060 


1.1245 


1.1343 


8 


0.83438 


0.8914 


1.0032 


1.0146 


10 


0.73731 


0.8082 


0.9134 


0.9259 



TABLE I: The values of ai^' for /3 = 1/4, 1, 2, oo corre- 
sponding to instability of flat phase state, (accuracy of cal- 
culation is ±0.0001) 



prime denotes a derivative with respect to a) 
P > 1/2. Using Eq. l(^gal) we have 



E+1 _ tan(-/?i(a)) 
1- E ^ tanh(/?a) 



< 1 



i.e., if 



(30) 



and, consequently, E < inside of the interval < a < 

We have found that for any P > 1/2 there is an interval 
< a < Oc where the flat phase state /ic = is stable. If 
P < 1/2, then E > and the flat phase state /ic = is 
unstable. Thus, /3 = 1/2 is a threshold value, for which 
E = and Oc = 0. 

In a similar way, one can show that the flat phase state 
/ic = TT, corresponding to Eq. (|23b|l . is stable inside the 
interval 0<a<ac, if/?<l/2. The appropriate system 
defining Oc is: 



tan(/3ac) — 

tanh(ac + Pj) ~ 

tan(ac Pj) = 

Pn/2 = 



- tanh(/3i), 

-tan(/3j+i), j = 1,3, . . 
tanh(/3j+i), j = 2,4,. 



(31) 



The values of Oc corresponding to the instability of the 
flat phase state are calculated using Eqs. (|29|l and H31|) 
for P — |;,l,2,oo, and summarized in Tab. One can 
check that the values in this table are in accordance with 
those obtained earlier by direct numerical simulationiS 
for P = 1/2,1. The result for iV = 2 and /3 = oo co- 
incides with the one calculated earlier analyticalljiii*2£. 
We stress here that our present results are obtained for 
arbitrary edge facet length b (arbitrary /?) and have much 
higher accuracy (can be calculated much faster). 

The plot ai^\p) for different N can be seen in Fig. |31 
As ^ 0, the af^ oo while for > 2 the a[^^ 
approaches some finite value. Also note the following 



natural property which can be seen in Fig. |21 

-,(2) 



{N) 



(1) 



(0). The fact that ac oo for /3 ^ means that 



— TT — LJJ approaches the limit of all-7r LJJ, where 
vortex solutions are unstable and flat phase state wins. 
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FIG. 3: The dependences of crossover distances ai^\ a^c'^ 
and ai^^ on j3. These curves can be obtained by solving either 
stability problem [Eqs. jSHJ for 6 > a/2 or Eqs. for b < 
a/2] or AFM solution existence problem [Eqs. 14611 for b > a/2 
and Eqs. lg7|l for 6 < a/2]. 



1. Behavior of instability point at large N 



Thus, the first Eq. of the system 1)32(1 now reads 

1 



at 
2 24 



(34) 



Finally, the expression for Oc is 



2 N 



(35) 



It is valid for N > 24(1/2 - /?). 

Second, for finite /3 > 1/2 the reasoning is similar, and 
we get 



^ 2 N 



(36) 



This approximation is valid for N ^ 24(/3 — 
1/2) max(l,/32). 

In the case P oo (infinitely long edge facets ), 
tanh(/3ac) — > 1 and our approximation does not work. 
The reason is that (3j does not vanishes with increase of 
N in this case. 

Alternative, more complex but also more exact, deriva- 
tion of asymptotic behavior of Oc for N oo and any 
finite as well as infinite (3 is presented in appendix 151 



Using the systems (|^ and its possible to de- 
termine a ground state of a very large LJJ, i.e., when 
N ^ oo. 

First, let us consider the case (3 < 1/2. As we saw 
above, according to the system (|31|l when 0, the 

values of f3j ~ Oc. We write Pj = acipj, where tpj ^ 1- 
In this case we expand the equations of the system (|31f) 
in Taylor series for small Oc and discard all the terms 
smaller than 0(a|?): 

r P ^ -1^1 + lai^pi 

I ~ - fa^V^f+i, J = 1,3,... 

[ V'JV/2 = -1/2. 

Substituting expressions for one by one starting with 
the last equation and discarding all the terms smaller 
than 0{a^) we result in the final expressions for ^pj: 

I 



B. Existence of AFM ordered solution 

In this section we apply the general formulas derived 
in the Sec. El to calculate the crossover distance for 
AFM state in the LJJ with equidistant distribution of 
TT-discontinuity points. 

If one chooses a domain of parameters where AFM or- 
dered SF chain is present, vary a and plot the solutions 
/x(a;), he may notice that the amplitude of spatial oscilla- 
tions of n{x) decreases as a decreases, see Fig. 0] It may 
happen, as we show below, that the amplitude of ^(x) 
vanishes at some a = which may be larger than zero. 

In this section we re-introduce the crossover distance 
as a distance at which the amplitude of oscillations 
of iJ.{x) vanishes, i.e., fJ.{x) — > fic for a flc + 0. In 
fact the limiting value of the phase /ic can have only 
three different values: 0, 7r/2, tt. To prove this, we refer 
to the Eq. © and write the set of expressions for the 
distance a = \x{fj,j + l)—x{fij)\ (j = 1, . . . , TV— 1) between 
discontinuity points 



dv 



y/Cj - 2(-l)Jcosj^ 



dv 



y/Cj - 2(-l)Jcosiy 



(37) 
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FIG. 4: Graphs of amplitude decreases with decrease of a: (a) b = a/2, A'' = 4; a = 3/2, 1, 1/2, = 7r/2; (b) fe = a, A'^ = 4; 
a = 3/2, 1.25, 1.2, fi^ = 0; (c) b = a/4, AT = 4; a = 3/2, 1.3, 1.27, /ic = tt 



j = 1, . . . , A^. For the left edge interval it reads: 

Ml 

I \ , (38) 
J V L'O - 2 cos J/ 

Due to the symmetry conditions (|22() we consider here 
only the left half of the junction. Note that in spite of 
index j, both sides of the expression H37(l do not depend 
on the facet number j, since we consider equidistantly 
distributed discontinuity points. Consider the limit a — > 
Uc + 0, which means 

/ij~Mc + e/ij, (39) 

where e ^ 0. The integration intervals are proportional 
to e. The denominator can be approximated using ly — 
fic + si^ as: 

\/—2€Ks sin /ic — e'^Kc cos fj,c, 

where Kg and Kc do not depend neither on fic nor on e. 
If /ic = or /ic = TT, the denominator in Eq. H37|l is oc e, 
which results in Oc > 0. Otherwise the denominator is 
oc ^/e, so that Oc — 0. 

For all other values of /Zc ^ and /ic 7^ tt, the de- 
nominator in Eq. H37|l is constant and integral vanishes, 
i.e., flc = 0. It is interesting that "all other values of 
lie' essentially mean /ic = 7r/2, see appendix 1X1 for de- 
tails. This also allows us to make a quick conclusion 
that Oc = for odd A^. Indeed, for odd A^, due to the 
symmetry conditions (|22b|l H(^n+i)/2 = i'/2 for any a. 



Therefore, /^(a;) fJ-c = 7r/2 when a decreases. This 
automatically means that Cc = 0. 

Summarizing our findings, one can get the following 
possible values of fic depending on Uc and (3: 

1. /ic = tt/2 for odd A^, Fig. 2fa); in this case Uc = 0. 

2. Hc^O for even N and P > 1/2, Fig.g^b); in this 
case flc > 0. 

3. = TT for even A^ and (3 < 1/2, Fig. gfc); in this 
case Oc > 0. 

Since, for odd A^ the = is already known, we 
calculate Uc only for the last two cases. 

Since for AFM state one has the symmetry (|22a|l . So, 
it is enough to consider j = 0, . . . A'^/2 with condition 

A*(aA'/2) = ^J'{aN/2-l)■ 

Even N, /3> 1/2, = 0. Using Eq. ^ with fic = 
the Eqs. ijH) and ^ can be approximated as follows: 

^2{-iy +e^j:,, (40) 
l,i^ « 6y^(-l)^+iS, = (41) 

where we defined Sj as 

S,=-2^(-l)^/i2_^2_ (42) 

2=1 

It follows from the Eqs. (221 and gDJ) that > for 
odd j and Sj < for even j. 

In the limit e ^ the Eqs. l|S7|) and become: 
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I. 



dv 



dv 



, j = l,...,iV/2. 



(43) 



diy 



(44) 



For the sake of simplicity we introduce new variables — jlj / jli, = 1, and 



\ 



1=2 



2(-ip+'(-E(-i)^4' + i-f )• 



(45) 



The integration of Eq. H43|l separately for odd j < N/2, even j < N/2 and the integration of Eq. H44|) give the 
system of equations for ai^^ and miscellaneous variables : 



(N) ■ C7 + I • ^7 11- ^ 

a). ' = TT — arcsm -4: arcsm -r^ , odd 7 < — 



AT 

even 7 < — 



(46a) 
(46b) 
(46c) 



i?uen N , (3 < 1/2, /ic = tt. Using Eq. 139|) with ^c — t^ and following the same procedure, we arrive to the following 
system of transcendental equations which define ac'- 



Pa 



(N) 



arcsm , 
2 VCo 



i^^-* — t: ~ arcsin ^i^i — arcsin -1^, for even j < N/2 



(47a) 
(47b) 



al^) = In 



, for odd j < N/2 



(47c) 



Now we consider infinitely long JJ. Similar to the finite 
LJJ, flc is not zero only for even N. For calculation of 
the crossover distance in this case we use equations (|46ll , 
with limit /3 — > 00 and /zq = 0. This modifies Eq. H46a() 
for j = 1: 

a(^) = -vr- arcsin^. (48) 
4 ^ 

The rest of equations from the system (I46II stay the same. 

To give an example, the values of ai^^ calculated for 
(3 = 1, 13 — 1/4 and (3 = 00 are presented in Tab. ITTl 
Tab, mil and Tab. II VI respectively and are in accordance 
with Tab.n This technique of solving numerically a sys- 
tem of transcendental equations is rather effective and 



can be used to obtain the plots ah. {(3) shown in Fig. 
Note, that these plots exactly coincide with the ones ob- 
tained in section IIII Al using stability analysis for flat 
phase state. 

The coincidence of the crossover distances obtained in 
two different ways implies that the transition between 
AFM state and flat phase state at a = Cc happens be- 
cause AFM solution just cease to exist. It was believed 
before^^. that transition takes place because one of the 
states has lower energy. Now it is also clear why in Ref.— 
the hysteresis around ac was never seen! Hysteresis usu- 
ally takes place when one has two stable solutions having 
different energies. 

We can draw possible states of the system as a pitch- 
fork bifurcation diagram shown schematically in Fig. |S1 
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iV 




Solution of ESI 


2 


1.4639 




Co = 0.4392 


4 


1.1772 




Co = 0.5628,6 = 1.1469 


6 


1.0060 


Co 


= 0.6451,6 = 1.1807,6 = 1-3141 



TABLE II: The values of ai for f3 = 1 (accuracy of calcu- 
lation is ±0.0001). 



iV 




Solution of EZIl 


2 


2.9277 




Co = 1.34429 


4 


1.2546 




Co = 1.0513, C2 = 1.3734 


6 


0.9833 


Co 


= 1.0310, C2 = 1.2352, Ca = 1.3233 



TABLE III: The values of a^' for /? = 1/4 (accuracy of 
calculation is ±0.0001). 

At small a < flc the flat state is the only solution and 
it is stable. At a = Oc the flat phase solution looses its 
stability, and it is unstable at a > Oc as indicated by the 
dotted line. In the same time, at a = Oc two new solu- 
tions appear. Both correspond to AFM ordered chain of 
semifluxons but with different sign. 




facet length, a 



FIG. 5: The sketch of the bifurcation diagram which show 
the transition from flat phase state to the state with AFM 
ordered semifluxon chain. 



iV 


JV 


Solution of Eq. EHl and EHl 


2 


7r/2 




4 


1.3063 


C2 = 1.2266 


6 


1.1343 


C2 = 1.3290, Ca = 1.6058 


8 


1.0146 


C2 = 1.3772, Cs = 1.7641, C4 = 1.9065 



TABLE IV: The values of aj, for /3 = 00 (accuracy of calcu- 
lation is ±0.0001) 



IV. CONCLUSIONS 

We have studied analytically the ground states in a 0- 
TT-LJJ with different number of facets of the length a ^ 
A J. We have shown that in the general case there is a 
crossover distance ai^"^ such that if the facet length a < 
ac^"* , the system is in the flat phase state {fi = const) and 

contains no magnetic flux. In contrast, if a > ai^\ the 
ground state consists of fractional vortices, each pinned 
at the phase discontinuity point. There may be more 
than one such state, especially for large N, but we focus 
our attention on the most natural one — AFM ordered 
chain of semifluxons. The system chooses between flat 
phase state and AFM ordered chain of semifluxons not 
because of the energy competition as it was suggested 
earlier—, but because there is only one stable solution 
for given a, as shown in the bifurcation diagram Fig. [S] 
for a < ai^'' AFM ordered semifluxon solution does not 
exist, while a flat phase state exists and is stable; for a > 
ai^^ flat phase solution /i = const exists but is unstable, 
so the state is AFM ordered semifluxon chain. We have 
calculated the crossover distances ai^^ and summarize 
our results as follows. 

• For odd N , — 0, semifluxons are always present. 

• For even N Uc > 0- The dependences of ai^', ai^'' 
and ai^^ on b are shown in Fig.|3| In particular, for 
b — a/2 the ai^' = and semifluxons are always 
present, for all other 6, Oc > 0. 

Our calculations of Qc agree with previous numerical and 
analytical resultsiitiSiSS, but cover also the cases of larger 
A'', arbitrary edge facets length b and have much higher 
accuracy. We also show that in many cases the size of the 
edge facets b can drastically affect the state of the whole 
system, especially when 5«a/2or6^0 and N ~ 2, as 
can be seen from Fig. |21 

We stress that wc derived the position of bifurcation 
point flc approaching it from both flat phase state (from 
the left in Fig. O and from the state with AFM ordered 
chain of semifluxons (from the right in Fig. |SJ| , and we 
got the same results. In the first approach, the system 
of equations [Eqs. (j^H) for b > a/2 or (EH for b < a/2], 
that describe the (in)stability of the flat phase state, is 
particularly easy to solve numerically and the reader is 
encouraged to do so for his/her favorite values of N and 
b just setting proper seed value for Oc. Nevertheless, our 
derivation of more complex Eqs. H4t)|) for b > a/2 and 
Eqs. (|47|l for 6 < a/2, which describes the disappearance 
of AFM ordered semifluxon chain and gives the same val- 
ues of flc, is not useless. This approach, although more 
complex, allows to find the existence region for more com- 
plex semifluxon states like ttii; which will be discussed 
elsewhere using the results obtained here. 

We have also found that the crossover distance ac cx 
1/V n for large N, see Eqs. (gSJ and ^ or Eqs. l|B9|) 
and ljB10|l of appendix ^ Having a flxed, the longer 
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O-TT-LJJ (larger N) favors configurations with semiflux- 
ons and, therefore, with magnetic flux, while shorter LJJ 
(smaller N) favors the state without flux. Instead, if we 
fix the total LJJ length, the LJJ with smaller a (large 
N) will favor flat phase state, while the LJJ with larger 
a (small N) will favor the state with semifluxons. 

In the future, it is quite interesting to consider the 
possibility to have less natural states like ttii- This will 
correspond to the additional branches on the bifurcation 
diagram and there will be clearly a minimum distance 
ac(TTii) > ac(TiTi) (ac(TiTi) is the one found here) for 
which such a state is stable. For a > ac(tTii) there will 
be energy competition between various states, e.g., be- 
tween titi and ttii- Next, in terms of studying classical 
and quantum tunneling between various states like ttii 
and titi J it is interesting to consider how Qc depends on 
the applied magnetic field and bias current. 
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APPENDIX A: THE ONLY "OTHER" VALUE OF 

He — 7r/2 

Now we show that if Uc ~ 0, then fic = 7r/2. Let us 
subtract expression (^J for two end values of /Lt in the 
ri-th odd interval: /i = /i„ and /x = fJ-n+i- For them 
w ±V2, are defined by the Eq. 
0. Having this, we receive: 



2am[(a;„+i --0/a„,TO„ 
2am[(a;„ - a;*)/a 
2Ve 



dn 



0(6) 



where dn(a;,m) — y 1 — msin^ am(a;, to) is the Jacobi 
elliptic function. There is only one term of the order -^/e 
in this expression. Thus this term equals zero. Using 
expression (jHJ for a:* [x* = Xn — a„ F(/Xn/2, to„)] with 
substitution to„ w 2 and d am{^,m)/d£, = dn(^,TO) we 
get: dn[F(A*c/2, 2), 2] = which has the only solution 
inside of the interval < ^ < tt: fic — 7r/2. 



APPENDIX B: BEHAVIOR OF INSTABILITY 
POINT AT LARGE A^: ASYMPTOTIC 
RELATION FOR ARBITRARY /? 

The analysis in this section is based on the fact, that 
{Pj+2 — <C 1 for iV ^ oo. This allows us to 

approximate the function Pj of discrete parameter, j, by 



the pair of continuous functions and derive the first order 
ordinary differential equation with boundary conditions 
for one of them. Solving it one gets the implicit relation 
between Qc and N. Without loss of generality we assume 
that N/4: is odd in this section. 

First, we consider the case /3 > 1/2. Let us intro- 
duce two functions, corresponding to odd and even inter- 
vals: A{n) = P2n-i, B{n) = /32„, n = 1, 2, . . . , N/A and 
rewrite the system (|29|l in the following form 



A(l) 

tan [flc -f- A{n)] 
tanh [flc -I- B{n)] 

A{NIA - 2) 



The index n = 1,2..., 7V/4 - 2 in Eq. |jB2l) and n = 
l,2...,iV/4-l in Eq. (lB3ll . 

Now we may write A{n + I) A{n) + A' {n) , A'{n) < 
A{n), where by definition 



- arctan [tanh(/3ac)] , (Bl) 

- tanh [B{n)] , (B2) 

- tan [A{n + 1)] , (B3) 

^W4) = -y. (B4) 



A{n + An) - A{n) 
An ■ 



A'(n) = lim 



Express B{n) in terms of A{n) using Eq. ljB2|l and expand 
the r.h.s. of the Eq. (|B3() in series with respect to A' 
keeping only linear term in A' . Then we have (we now 
write "=" instead of "w"): 



A' = f{A) = 



cos A [sin flc — cos(ac + 2A) tanh Oc 



(B5) 



cos(ac + A) — sin(ac -|- A) tanhoc 
A{1) = - arctan [tanh(/3ac)] , A{N/4) ^ -ac/2.(B6) 

where A' is positive, since A{n + 1) > A{n). The later is 
consequence of the Eqs. (|B2|I . (|B3|I and the property of 
the increasing concaved down function /: f{xi -\- X2) < 
f(xi) -\- f{x2). Thus the above equation may be formally 
integrated: 



-ae/2 



A{1) 



d(3 

W) 



(B7) 



Although integration can be done explicitly, we stay with 
this symbolic form for the sake of simplicity. 

Similarly, for /3 < 1/2 we get the same Ea. (|B7|) with 
different function / and boundary condition A{\) = 
— atanh [tan(/3ac)]: 



I{A) 



cosh A [sinhflc — cosh(ac -t- 2 A) tanoc 
cosh(ac + A) + sinh(ac + A) tanoc 



(B8) 



with negative A' . 

Ea. l|B7p gives us essentially function N{ac) rather then 
desirable ac{N). Fortunately N{ac) may be simply in- 
verted since Oc «C 1 everywhere in these calculations. 
In fact, one can expand l.h.s. of ljB7|) in powers of Oc, 
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10 100 

number of facets, N 



FIG. 6: Asymptotic behavior of ac{N) for /? = 1. Solid bold 
line is received using Eq. IIB7t ; dotted line corresponds to the 
Eq. IIBQfl : solid thin line is result of numerical solution of the 
system 11291 . 



In particular, for /3 — s- cx) we have g = n\/3/2. 

Thus for large N and any length 6 = /3a of the end 
facets we have derived asymptotic relation Oc ~ N~^/'^, 
which is in agreement with equations of the Scc lIII ATI 

Now we derive applicability condition for the equations 
of this section. Eq. (|B7p can be used for such N that 



max A' 



A 



< 1. 



(Bll) 



Investigation of the series in Cc of the ratio (max^')/yl 
shows that it does not have local extremum. Thus its 
maximum value can be taken only at the end points of 
the integration interval: A = — ac/2 or A = A{1). One 
can show that this point is A{1) for all finite values of /3 
as well as for infinite (3. 

Using leading (linear) term of the series of A' in we 
get from Eqs. (|B5|I and 



keeping only the leading term, which is of the order ^ . 
Thus 



A' 
A' 



-2acSm^{A), /3 > 1/2, (B12) 
-2acSinh2(yl), /3 < 1/2. (B13) 



where 



N '■ 



(B9) 



g{l3) = V3arctan ^3(2/3-1) , P > ^, (BlOa) 
g(/3) = \/3arctan \/3(l - 2/3) , ^ < i. (BlOb) 



Together with relation (|B9|I and the fact that A{1) « 
-/3ac for Pac < 1, -7r/4 < A{1) < for (3ac > 1 and 
A{1) = — 7r/4 for /3 oo, we arrive to 

N > 8g, Pac < 1, (B14) 
N > 2V2g, /3 ^ oo or /3ac > 1. (B15) 
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